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Plasma 


What  is  Plasma? 

a  A  plasma  is  a  gas  in  which  an  important  fraction  of  the  atoms 
is  ionized,  so  that  the  electrons  and  ions  are  separately  free. 
How  to  model  Plasma? 

a  We  need  to  keep  track  of  its  evolution  both  in  spatial  and 
velocity  space. 

Modeling  methodology: 
a  Fluid  methods 
a  Kinetic  methods 
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Driven  Applications 


o  Electric  Propulsion 
9  Laser-Plasma  Interactions 
9  Hypersonic  Re-entry 
9  Material  Processing 
9  Astrophysics 
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Governing  Equations 


Starting  from  the  Boltzmann  equation  of  the  distribution  function 
of  a  multi-component  plasma 

df 

-q—  +  v  '  V xfa  +  3  ■  V vfa  =  Ca[) (fa)  ( 1 ) 

P 

Complete  solution  requires  solving  the  convection  equations  for 
each  plasma  component  in  6  dimensional  coupled  with  the  source 
terms  on  the  RHS  —>  highly  computationally  and  memory 
intensive,  model  reduction  required  — >  fluid  model 
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Fluid  Description 


Assumption 

a  Velocity  distribution  function  is  close  to  a  Maxwellian 
distribution. 

Approach 

0  Model  the  plasma  using  the  fluid  dynamical  equations  (Euler, 
Navier-Stokes) 

a  Thermochemical  non-equilibrium  environment  can  be 
characterized  via  multiple-temperature  model 

0  Eletromagetic  field  effect  can  be  modeled  using  the  MHD 
approximation. 
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Fluid  Equations 


Two-temperature  plasma  (T/,  —  Te ) 
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Solution  method: 

a  Finite  Volume  method  for  hyperbolic  conservation  laws 
9  Source  terms  are  solved  by  using  operator  splitting  technique 
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Kinetic  Description 


Assumption 

a  High  temperature 
a  Low  density 
Approach 

a  Treat  the  plasma  as  a  set  of  pseudo-particles  moving  in  the 
electromagnetic  field  computed  on  a  grid  (Particle  |n  Cell). 

a  The  plasma  is  then  modeled  by  intergrating  the  equation  of 
motion  (Lorentz  equations)  for  each  particle  inside  the  field. 
The  field  is  also  updated  by  solving  the  Maxwell  equations. 

a  Ionization  can  be  modeled  by  particle  collision  (e.g. 
Monte-Carlo  collision) 
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Particle-in-Cel I  Method 


Particle  equation  of  motion 


q(£  +  v  x  B) 


(4) 


Field  Equation 

o  Electrostatic:  Poisson  equation 

V20  =  —  (Z;/7;  -  ne)  (5) 

eo 


where  £n  =  — Vn0 

9  Electromagnetic:  Maxwell  equations 
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Current  Focus 


Current  Focus 

o  Basic  framework  for  Fluid  (CFD)  and  Particle  (PIC) 
simulations  on  GPU 

o  Fluid:  CFD  solver  with  reaction  kinetics 
»  Particle:  Electrostatic  PIC 
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Graphic  Processing  Unit 


What  is  GPU? 

a  Graphic  processing  units  containing  a  massive  amount  of 
processing  cores 

o  Designed  specifically  for  graphic  rendering  which  is  a  highly 
parallel  process 

Why  GPU? 

9  GPU  is  faster  than  CPU  on  SIMD  execution  model 
9  GPU  is  now  very  easy  to  program 
9  GPU  is  much  cheaper  than  CPU 
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GPU  versus  CPU 


2003  2004  2005  2006  2007  2008  2009  2010 


Figure:  Single  and  double  precision  floating  point  operation  capability  of 
GPU  and  CPU  from  2003-2010 
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GPU  Programming 


Programming  languages  for  GPU:  CUDA,  OpenCL, 
DirectCompute,  BrookGPU,  ... 

CUDA  is  the  most  mature  programing  environment  for  GPU. 
o  similar  to  C/C++ 
o  support  00  features 
9  easy  to  debug 
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GPU  Programming  Model 


9  Based  on  grid  and  thread  blocks 
9  Execution  instruction  called  kernel 


Multithreaded  CUD*  Program 

WBi  MB  MB  MB 
m  mb  mb  mi 
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GPU  Implementation 


CFD: 

9  Map  the  computational  domain  with  CUDA  grid 
o  All  the  computational  cells  can  be  processed  in  parallel 
PIC: 

9  Associate  particles  with  threads  and  have  each  thread 
integrate  the  equation  of  motion. 

9  Associate  field  quantities  of  the  cells  with  threads  and  have 
each  thread  update  the  field  value  of  the  cell. 
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CFD 


Forward  step  -  Mach  3  flow  over  a  step 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Plasma  Simulations  on  GPUs 


Introduction 
Plasma  Simulation 
GPU  Implementation 
Results 

Conclusion  and  Future  Works 


CFD 


Backward  step  -  Mach  2.4  shock  diffraction 
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CFD 


Rayleigh-Taylor  Instability  -  Acceleration  of  a  heavy  fluid  to  a 
lighter  fluid 


GPUs 
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PIC 
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Performance 


Without  chemical  kinetics 
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Performance 


With  chemical  kinetics 
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-  9  species,  38  reactions 

-  36  species,  308  reactions 
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Accomplishment: 

o  Basic  CFD  framework  for  fluid  simulation  and  PIC  framework 
for  electrostatic  field  on  the  GPU 

9  Performance  obtained  in  both  cases  are  very  promising 
Future  Works: 

9  Plasma  ionization  kinetics  (CR  model)  for  fluid  modeling. 

9  Monte-Carlo  collision  for  PIC  (PIC-MCC) 

9  Multi-fluid  and  hybrid  modeling 

a  Model  each  component  of  the  plasma  by  different  methods 
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